Perturbation Theory for Fractional Brownian Motion in Presence of Absorbing 

Boundaries 



Kay Jorg Wiese 1 , Satya N. Majumdar 2 , and Alberto Rosso 2 

1 CNRS-Laboratoire de Physique Theorique de VEcole Normale Superieure, 24 rue Lhomond, 75005 Paris, France 

2 CNRS-Laboratoire de Physique Theorique et Modeles Statistiques , Universite Paris-Sud, 91405 Orsay, France. 

Fractional Brownian motion is a Gaussian process x(t) with zero mean and two-time correlations 
(x(ti)x(t 2 )) = D (ti H + t\ H - \t\ - t 2 \ 2H ), where H, with < H < 1 is called the Hurst exponent. 
For H — 1/2, x(t) is a Brownian motion, while for H ^ 1/2, x(t) is a non-Mar kovian process. 
Here we study x(t) in presence of an absorbing boundary at the origin and focus on the probability 
density P+(x, t) for the process to arrive at x at time t, starting near the origin at time 0, given that 
it has never crossed the origin. It has a scaling form P+(x,t) ~ t~ H R+{x/t H ). Our objective is to 
compute the scaling function R+(y), which up to now was only known for the Markov case H — 1/2. 
We develop a systematic perturbation theory around this limit, setting H = 1/2 + e, to calculate the 
scaling function R+(y) to first order in e. We find that R+(y) behaves as R+(y) ~ as y —> (near 
the absorbing boundary), while R+(y) ~ y 1 exp(— y 2 /2) as y — >• oo, with (j) = 1 - 4e + 0(e 2 ) and 
7=1 — 2e + 0(e 2 ). Our e-expansion result confirms the scaling relation = (1 — H)/H proposed in 
Ref. [29]. We verify our findings via numerical simulations for H = 2/3. The tools developed here 
are versatile, powerful, and adaptable to different situations. 



I. INTRODUCTION 

Survival of a species of bacteria, translocation of DNA 
through a nano-pore, and diffusion in presence of an ab- 
sorbing boundary are only few out of many situations, 
where the central question is the survival, or persistence 
of the underlying stochastic process. More precisely, per- 
sistence, or survival probability S(t) of a process is the 
probability that the process, starting from an initial pos- 
itive position, stays positive over a time interval [0,t]. 
For many stochastic processes arising in non-equilibrium 
systems, persistence decays as a power law S(t) ~ t _6> , 
where is called the persistence exponent [I] . For a sim- 
ple Markov process such as one-dimensional Brownian 
motion, = 1/2 [2 . On the other hand, the exponent 
is non-trivial whenever the process is non-Mar kovian, 
i.e., has a memory. In addition to theoretical studies (for 
a brief review see [3 ), the exponent has been mea- 
sured in a number of experiments [4-10 . Even for Gaus- 
sian non-Markovian processes, is non-trivial [11]. For 
the latter processes that are close to a Markov process 
(i.e., whose correlators are close to that of a Gaussian 
Markov process) the exponent 6 was computed perturba- 
tively [12] [13] . This perturbation theory has been used 
for various out-of-equilibrium systems, as the global per- 
sistence at the critical point of the Ising model in d = 4— e 
dimensions [H] , in simple diffusion close to dimension 
[15], and in fluctuating fields such as interfaces [T6UT8] . 

A quantity that contains more spatial information than 
persistence S(t) is the probability density P+(x,t) of the 
particle at position x and at time t, given that it has 
survived (stayed positive) up to time t. To investigate 
P+ one can equivalently think of a process on the 
positive semi-infinite line [0, oo] with absorbing boundary 
condition at the origin x = (see Fig. [TJ. The question 
is, how does P + {x,t) depend on xl In other words, how 
does the presence of an absorbing boundary at the origin 
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FIG. 1: (Color online) The fractional Brownian motion dis- 
cussed in the main text. 



change the spatial dependence of the probability density 
of the particle at time tl In particular, it is clear that 
P+(x, t) must vanish as x — » and x — >• oo. But how do 
they vanish there? One of the main messages of our pa- 
per is that for generic non-Markovian processes, P+(x,i) 
vanishes near its boundaries at x = and x — » oo in a 
non-trivial way, characterized by non-trivial exponents. 

As the persistence S(t), the probability P+(x, t) can be 
computed exactly for a Gaussian Markov process, as e.g. 
one-dimensional Brownian motion. For non-Markovian 
processes, even if they are Gaussian, P + (x,t) was not 
known. In this work, we consider P+(x,t) for a class of 
one-dimensional Gaussian processes known as fractional 
Brownian motion (fBm), which are parametrized by their 
Hurst exponent H, with < H < 1. The case H = 1/2 
corresponds to ordinary Brownian motion, which is a 
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FIG. 2: Left: Translocation of a polymer chain through a 
pore. Right: The translocation coordinate s(t) denotes the 
number of the monomer that is crossing the pore at time t. 



Markov process, while for H ^ 1/2 the process is non- 
Markovian. The purpose of this paper is to develop a 
systematic perturbation theory to compute P+(x,t) for 
non-Markovian fBm's with H = 1/2+e, where e is the ex- 
pansion parameter for the perturbation theory. Here we 
present the result for P+(x,t) to 0(e). It can be written 
as a combination of special functions, i.e. error and hy- 
pergeometric functions, see Eq. ( flQ] ) . To our knowledge, 
this is the first systematic (exact up to 0(e 2 )) calculation 
of P+(x, t) for fractional Brownian motion with H ^ 1/2. 

Before detailing our results, let us position them into 
a broader context: Fractional Brownian motion with 
H 7^ 1/2 is relevant for polymer translocation through 
a nanopore. Consider a polymer chain composed of TV 
monomers passing through a pore (translocation) from 
left to right, as drawn on Fig. |2j The dynamics of this 
translocation process has been investigated intensively 
due to its central role in understanding, e.g., viral in- 
jection of DNA into a host, or RNA transport through 
nano-pores, and mastering such applications as fast DNA 
or RNA sequencing through engineered channels [TM22] . 
The translocation coordinate s(t), namely the label of the 
monomer crossing the pore at time t, is key to quantita- 
tively describing the translocation process [23- 26 , which 
begins when s = 1, and ends when s = TV, i.e., when the 
first and the last monomer of the chain enter the pore, 
respectively, see Fig. [2] For large A/", when the transloca- 
tion is not yet complete, one can view s(t) as a stochastic 
process on the semi-infinite line with absorbing boundary 
conditions at s = 0. The absorbing boundary at s = 
models that if the chain falls back to the left, i.e., on 
the starting side, it will diffuse away and not try again. 
The quantity P+(s(t) = x,i) then represents the prob- 
ability that x monomers have translocated to the right 
at time t. To model the process s(t), one observes the 
following facts: (i) scaling arguments and numerical sim- 
ulations show that s(t) is subdiffusive [27]; (ii) in absence 
of boundaries, numerical simulations indicate that s(t) is 
a Gaussian process [28]. Based on these observations it 



was proposed in Ref. [29 that a good candidate for s(t) 
is a fractional Brownian motion with H = 1/(1 + 2z/), 
where the exponent v describes the growth of the radius 
of gyration with the number of monomers (R g ~ N u ) 
[30] . Thus for ^^1/2, 1/2 and hence s(t) is generi- 
cally a non-Markovian process, with absorbing boundary 
conditions at s = and at s = N. Here we consider 
the limit of N — >• oo. Thus our results for P+(x,t) of 
a fBm with H ^ 1/2 are directly relevant for polymer 
translocation. 

Directions for further applications are numerous: Re- 
cently a relation was established between the statistics 
of avalanches associated with the motion of a driven par- 
ticle in a disordered potential and persistence proper- 
ties of the latter [31] . Higher-dimensional generalizations 
are avalanches of extended elastic objects, for which sys- 
tematic field-theoretic treatments exist [32H35] . In few 
cases, no-hitting probabilities can be calculated for ex- 
tended (non-directed) objects, as self-avoiding random 
walks avoiding extended objects [36]. Other approaches 
use real-space renormalization [371 13H] • 

This article is organized as follows: Since some of the 
computations are rather technical, we first provide in Sec- 
tion [IT] a brief summary of the main definitions and our 
principal results. In Section [Hi] we introduce basic nota- 



tions and reproduce the known results for H = £ . Section 
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IV| explains the basic ideas of our perturbative approach, 
sketches the calculation, and discusses some of the sub- 
tle points. Our predictions are compared to numerical 
simulations in Section (Vl Conclusions are presented in 



Section [VT} More technical points are relegated to two 
appendices: In Appendix [A] the correction to the action 
is derived. Appendix [B] contains the explicit calculation 
of the perturbation theory. Finally, Appendix [C] reviews 
the arguments for the scaling law <p = (1 — H)/H. 



II. SUMMARY OF DEFINITIONS AND MAIN 
FINDINGS 



Consider a particle, located at time t = at the origin 
x = and free to propagate on the real axis. For Gaus- 
sian processes, the probability to find the particle inside 
the interval (x, x + dx) at time t is given by 



P(x, t) dx 



1 



e 2 <- 2 (*)> dx , (1) 



where (x 2 (t)) is the particle's mean square displacement. 
A natural scaling variable is 



(2) 



and most of the properties of the process are a function of 
this single variable. For example, the distribution prob- 
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ability in Eq. becomes 
P(x, t) dx 

R(y) = - 



= R(y)dy 
1 _Z 
/2tt 



(3) 
(4) 



In many problems the motion is confined to an interval, 
finite or semi-infinite. In presence of absorbing bound- 
aries, the probability distribution of the particle position, 
subject to the condition that the particle has survived, 
has no longer a simple Gaussian form since it has to van- 
ish at the boundaries. However, one can still express it as 
a function of the sole scaling variable y defined in Eq. ([2| , 
where (x 2 (t)) is the particle's mean square displacement 
in the unconstrained (without boundaries) process over 
the full real line. In particular, here we discuss the case 
where the particle can move on the positive semi-axis 
and is absorbed whenever x(t) < 0. We call P + (x,t) and 
R+ (y) with y given in Eq. Q the normalized probability 
distribution and the scaling function of the problem in 
presence of an absorbing boundary at the origin, 



P+(x,i) dx = R+(y) dy 



(5) 



In contrast to the free case, the functional form of R+(y) 
is not the same for all Gaussian processes, but depends 
on the precise nature of the latter. Here we study a par- 
ticular class of processes, the fractional Brownian motion 
(fBm), for which the autocorrelation function in absence 
of boundaries is 



(x(ti)x(t 2 )) =D(4 



2H 



t 



2H 



t 2 



\2H\ 



(6) 



where H with < H < 1 is the Hurst exponent. For 
H = 1/2, the fBm identifies with Brownian motion 



(x(h)x(t 2 )) = 2£>min(ti,t 2 ) 



(7) 



where D is the diffusion constant. Note that only for 
H = 1/2, the Gaussian process x(t) is Markovian. For 
other values of if, the process is non- Markovian. 

For Brownian motion (if = 1/2), the form of R+{y) 
can be obtained using the method of images (see Section 



Rf{y) = ye 



(8) 



The superscript (0) identifies the case if = 1/2. For other 
values of if, due to the non-Markovian nature of the 
process, the method of images no longer works and the 
computation of R+(y) becomes a challenging problem. In 
this paper we compute this function, using a perturbative 
approach for if = 1/2 + e, to first order in e. The final 
result is 



R+(y) 

W(y) 



- »(°) 



eW(y) + 0(e 2 
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t/ 2 )erfi(-M 



(9) 
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7E ] - 3y 2 , (10) 




FIG. 3: (Color online) The function W(y) denned in Eq. (10 
The solid blue line is the result of the series expansion (11 
while th e da shed red line is the result of the asymptotic ex- 
pansion O. We note that 1^(1.01694) = 0. 



where 7e is Euler's constant, 2^2 a hypergeometric func- 
tion, and erfi the imaginary error-function. We can write 
a convergent series-expansion 



W(y) = % 4 J2 

00 

-E 



2 n n! y 2n 
(2n + 4)! 

^ 2 j-ryn+i 
(2n- l)(2n + l)n! 



+ (^ 2 -2) [ln(2^ 2 )+ 7E ] -3y 2 (11) 

where each line is equivalent to the corresponding line in 
Eq. (10). Both sums converge for all y, but problems of 
numerical precision appear for y > 7. In that region, one 
can use the asymptotic expansion 

W(y) = l- 7E -ln(2 y 2 ) + ^-^ + A +0(y -8) (12) 



At y = 7, the difference between (11) and (12) is smaller 
than 10~ 6 . 

We obtain, at first order in e, the asymptotic expan- 
sions of R+(y), 

R + {y) ^ </[l-4elny 2e( 7E + ln2) + ...] 



R+(y) ^^ye-y ' 2 [1 - 2ehiy + e(l - In 2 - 7 e)] 



(13) 



These asymptotics can be recast into 

R+{y)~y* for y^0 

R+(y) ~ y 1 e 2 for y — >• 00 , 



(14) 



where the two exponents (j) and 7 are at first order in e 
given by 

= l-4e + 0(e 2 ) , 7 = l-2e + 0(e 2 ) . (15) 
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In a recent publication [29] (reviewed in Appendix [C| , 
a general scaling relation, valid for arbitrary self-affine 
processes with stationary increments, was proposed be- 
tween the exponent </>, the persistence exponent and 
the Hurst exponent H, 



r H ' 

For £Bm, it is known rigorously that ■ 
result predicts that for fBm, 

. 1 - H 



H 



(16) 

l-H [16 . This 
(17) 



One of the objectives of this paper was to verify this 
scaling relation up to 0(e) in a perturbation theory 
around H = 1/2. Using H = 1/2 + e, one expects 
(j) = (1 - H)/H = 1 - 4e + 0(e 2 ) for fBm. This is in 



agreement with our result (15), putting the scaling argu- 
ments on a firmer footing. 

It is interesting to note that the scaling function R+ (y) 
given in Eq. Q has, at least to 0(e), the same lead- 

_ 2 ley 

ing large-?/ behavior ~ e~ y 1 as in the unconstrained 
case Q. This behavior can be understood by a simple 
heuristic argument: far from the boundary the process 
is not "aware" of the latter. Our calculation reveals that 
the process nevertheless knows about the boundary, and 
R+(y) has a subleading power-law prefactor y 1 where 7 
is a new (independent) exponent, whose result to order e 
is given in Eq. ( 15 ). 



Our analytical results are then verified via numerical 
simulations for H = 2/3. 



III. PRELIMINARIES: BROWNIAN CASE 

(H = 1/2) 

To simplify notations, we set D = 1 in the following. 
The final result expressed in the variable y, is of 
course independent of this choice. 

The spreading of a Brownian particle is given by the 
Fokker-Planck equation 

d t Z^\x ,x,t) = d 2 x Z^(x ,x,t) (18) 

Z^\x ,x,t = 0) = S(x-x ) (19) 

The propagator Z+\xo,x,i) times dx gives the proba- 
bility to find the Brownian particle inside the interval 
(x, x -f dx) at time t, knowing that the particle was at xq 
at time t = 0. With absorbing boundary conditions at 
the origin we have, using the method of images 

1 



zf(x^x,t) = 



-{x-xo) 2 /At 



-(x+xoY/At 



(20) 

This propagator is not a probability distribution because 
it is not normalized. Its normalization, the so-called sur- 
vival probability, 



S(x ,t) 



r>00 

I &x Z+\xQ,x,t) 
Jo 



(21) 



gives the probability that the particle is not yet absorbed 
by the boundary at x = 0. The survival probability van- 
ishes when xq — >• 0; however, in that limit, the probability 
distribution for the non-absorbed particles remains well- 
defined: 



Pi°\x,t) = lim 



Another quantity with a finite limit for xq = is 



Z ( +\x,t)= lim —Z y +\xo,x,t) = 

x ^0 Xq 



(0), 



xe 4t 
2v^ 3 / 2 



This allows to write the probability P+\x,i) as 



ifW) 



Zf\x,t) 



2t 



(22) 



(23) 



(24) 



J °° dx Z^ (x, t) 

Using in Eq. ( |24| ) the scaling variable defined in ([2|, 
y = x/y/2t, we recover is\. Eq. (24) is simpler than 



Eq. (|22j) because the xq dependence is discarded from 
the beginning. We will use this definition to compute 
Z+(x,t) for H = 1/2 -he. 



IV. PERTURBATION THEORY (H / 1/2) 

The process x(t) is Gaussian for all values of H, but 
it is Markovian only for H = 1/2. For all other val- 
ues of H, the process is non-Markovian and this makes 
the problem difficult to solve. Our idea is to expand 
around H = 1/2. In a first step, we construct an action, 
which calculates expectation values of the Gaussian pro- 
cess x(t), with bulk expectation values (|6|. In a second 
step, we obtain the propagator with absorbing boundary 
conditions at x = 0. In a third step we calculate the 
probability P+(x,t) perturbatively, using the action con- 
structed in step 1. In the fourth step, we put together 
all pieces and interpret our result. 



A. Step 1: The Action 

For all H, x(t) is a Gaussian process, therefore the 
statistical weight of a path x(t') without any boundary 
is proportional to exp(— S[x\) where the action S[x] is 
quadratic in x and given by 



S[x] 



= I dh [ 
Jo Jo 



dt 2 ^x(t 1 )G(t 1 ,t 2 )x(t 2 ). (25) 



Note that we use standard field-theoretic notation, noting 
f(x) a function of the variable x, and S[x] a functional, 
depending on the function x(t / ), with < t r < t. 

The kernel C?(ti, t 2 ) of the action is related to the auto- 
correlation function of the process via 



G-\t u t 2 ) = (x(h)x(t 2 )). 



(26) 
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For H = 1/2, the action is simple. In this case, setting 
D = l, 

[G(°>]- 1 (*i,*2) = (x(h)x(t 2 )) = 2min(ti,t 2 ). (27) 



Using the result (|A7j) in Eq. (|25|, we recover the standard 
Brownian action 



S^[x) = - A j\t' {d t ,xf. 



(28) 



For a generic value of H the kernel G{ti,t 2 ) becomes 
non-local. For H = 1/2 + e one can write 



S[x\ =S (0) [x}+eS (1) [x] 



(29) 



where [x] is the action ( |28| and [x] has been com- 
puted in Appendix [A] 



o Jt 



dti / dt 2 



d tl x(ti)dt 2 x(t 2 ) 



\h-t 2 \ 



-2«S (0) [a;](l+lnr) . 



(30) 



Note that we have introduced a regularization for coin- 
ciding times t\ = t 2 — >• In \ ti — t 2 \ = lnr where r > is 
the UV cutoff. A first-principle definition would necessi- 
tate a discretization in time. It is however sufficient to 
check that the law ([6| is correctly reproduced, and that 
the final result is cutoff independent. 



B. Step 2: The Propagator with an Absorbing 
Boundary 

For a generic value of H, the propagator Z+(xo,x,t), 
denoting the probability that the particle reaches x at 
time t, starting from xq at time 0, and staying positive 
over the interval [0, t], can be written using standard path 
integral notation as 



Z+(x ,x,i) 



L 



x{t)—x 



x(0)=x o 



V[x]e~ s ^e[x] . (31) 



Here Q[x] is an indicator function that is 1 if the path 
x(t') stays positive over the interval [0, t] and other- 
wise. The action S[x] is given in (25). In the limit 
xq — 0, we expect, as in the Brownian case (H = 1/2), 
the propagator to vanish as Xq° where the yet unknown 
exponent </>o depends on H. Note that for H = 1/2, 
0o = 1 (see Eq. (|23|)). For H = 1/2 + e, we expect that 
cj)Q = 1+ai e+0(e 7 )7 where ai is yet unknown. Analogous 
to Eq. S for H = 1/2 we define Z+(x,t) as 



Z+(x,t)= lim 4r / ° ^We-^eb] . 



(32) 



Using the expansion of the action given in Eq. (29) and 
0o = 1 + clic, we write to leading order in e 



Z+(x,t) 
= lim - 

^o^0 x. 



1 



1+aie 




a:(0)=a:o 



V[x] (l-e<S (1 )[x])e- 5(0) WeH 



lim { zi°^(x, t) [1 — aieln(xo)] 



£c(t)=£C 



P[:r]«S (1) [ 



_ 5 (0) 



x Jx(0)=x 

zf\x,t) + eZ^\x,t) , 
AO), 



[x] e[x]} 



(33) 



where Z+ } (x,t) is defined in Eq. (23) and Z+\x,i) is 

-i px(t)=X 

— / V[x]S {1) [x]e 

x J x (0)=xo 



si0) ^G[x] 
-a 1 \n(x )Z ( ^\x,t)} (34) 

We will see that for O = 1 — 4e, i. e. ai = —4, Z^(x,t) 
is independent of xq. 



C. Step 3: Calculation of Z^\x,t) 

The main achievement of this paper is the calcula- 
tion of Z+\x,i) defined in Eq. (34). This calculation 



is rather involved, both conceptually and technically. 
Therefore, we will relegate several technical calculations 
to Appendixpl Eq. p4|) can be devided into three pieces: 



Z%\x,t) = Z*(x,t) 



(35) 



- lim 

x ^0 



Z+(xq, x, t) — a\ hi(xo)Z+\x, t) 



Z*(x,t) = 2(1 + lnr) (36) 

-i px(t)=X 

x lim — / V[x]S^[x]e- si0) ^e[x] 

x ^0 x J X (0)= XO 

rt pt 

yB< ~ ~ * - 1 1 1 dhdt 2 (37) 



Zf(x ,x,t) 



JO 



J X 



x — / V[xY—— — — — e 

XQ J x (0)= Xo \tl - t 2 \ 



e[x] 



The first term, Z+(x,t) is simple, and is evaluated in 
appendix B 1 We now come to the evaluation of the 
contribution Zf(xo,x,i), defined in Eq. (37). In Fig. [i] 
we show a path which contributes to Zf(xo,x,t). The 
sum of all these paths is a product of transition proba- 
bilities. Explicitly, it reads, ordering t\ <t 2l which gives 
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space 




* time 



FIG. 4: Graphical repr esen tation of the path-integral for 
Z+(x ,x,t) given in Eq. (34). 



an extra factor of 2 compared to (37): 



Z+(xo,x,t) 



/ / / 

X!>0 xt>0 x 2 >0 x 2 >0 

Z^\x 0j x u t 1 )D(x u x 1 )- + ( Jl - l2 ' f2 tl) 



1*2 "tl 



xD(x 2 ,x 2 )Z+\x 2 ,x,t - t 2 ) 



(38) 



Z^\xi,x ,ti), z\_\x 1 ,x 2 ,t 2 -t 1 ) and z\_\x 2 ,x,t-t 2 ) 
are defined in (20). The factors D(x\,xi) and D(x 2l x 2 ) 
take into account the terms d^x^ti) and dt 2 x(t 2 ) in the 
action Si [x] . 



D(xi,xi) = lim - Z^ (x\, x\, dt) 
v ' J dt^o dt + v 7 

, _ v _ (*i-*i) 2 
(xi — xi) e 4d* 



lim , . 

dt^o dt V27rdt 

-2^(xi -xi). 



(39) 



Finally, we have set t\ = t\, and t 2 = t 2l since we 
have taken the limit of their differences to 0. In or- 



der to perform the six integrations in Eq. (38) it turns 
out to be convenient to evaluate its Laplace transform, 
Z?(xq,x,s). From now on, we will always denote with 
f(s) the Laplace transform of a function /(£), defined as 



m := / 
Jo 



dtt 



■fit). 



(40) 



This Laplace transform leads to two important simplifi- 
cations: The first simplification is that now the nested 
time-integrals over t\ and t 2 become a product. To see 
this, we remind that if f\ and f 2 are two functions which 
depend on t, then the Laplace transform of their convo- 



lution is simply the product of their Laplace transforms, 



die 



■*2) 



/ d*i/i(*i)/ 2 (*-*i) 
Jo 

pOO pOO pOO 

dt dti / dt 2 5(t-ti 
Jo Jo Jo 

xfi(h)e- tlS f 2 (t2)e- t2S 

- ncx) "1 r poc 

/ dti/i(ti)e-* ia / dt 2 f 2 (t 2 )e- t2S 
Jo J Uo 

fi(s)h(s) . 



(41) 



This consideration generalizes to 3 and more times. 
We obtain for the Laplace transform of (38) 



Z^(x ,x,s) 



X 



(x ,x u s)Z^\x 2 ,x,s) 

(0)/ 



xd Xl d x 



xi>0 x 2 >0 

Z^>(x u x 2 ,t) 



die" 



t 



(42) 



The second simplification is even more important, and 
is most easily understood on the example of the bulk 
propagator 



Z (0 \x,y,t) := 



Its Laplace-transform is 



-(x-yf/4t 



Z^\x^s) = ^ r e-^-y\ . 
2v s 



(43) 



(44) 



While integrals over x > involving (B19) give error- 
functions, which are hard to integrate further, the same 



integrals over (44) remain similar exponential functions; 
the only complication is that one has to distinguish be- 
tween x smaller or larger than y. 

To evaluate (42 ), we now have to calculate the Laplace- 



transforms of its factors: 



p -y/s\x-y\ _ p -Vs{x+y) 

ZV {x,y,s) = ^= . (45) 



Finally, the term in brackets in Eq. ( [42] ) can be rewritten, 
using a Fourier decomposition for Z+ (x 2: xi : t), as 



dte 



t 



dk f°° , e 

— / dt 

2Wo t 



(s+k 2 )t 



ik(x!-x 2 ) _ e ik(x!+x 2 ) 



dk 

2^ 



e %k{x 1 -x 2 ) _ e ik(x 1 +x 2 ) 



[ln([ S +fc 2 ]r)+ 7E ] 



(46) 



Note that the time integral in the second line of Eq. (46) 



is diverging at small times. Since the path integral is 
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defined as discretized in time, a natural approach con- 
sist in discretizing this integral, with a step-size r. This 
would indeed be the only possible approach for stronger 
divergences, like 1/t 2 . However, since our integral is only 
logarithmically diverging, we can take an easier path, by 
using a small-time cutoff r: 



I 



oo e -( s +k 2 )t 



dt 



oo e _( a+ ^)t 



(47) 



t 



-dt 



■ln([5 + /c 2 ]r)-7E + 0(r) . 



We note that the regularization by discretization gives 
the same result apart from the term — 7e- We will check 
later that it only contributes to the normalization, which 
will drop from the final result. 



Collecting the results of Eqs. (45) and (46) in Eq. (42), 



and doing the remaining space-derivatives we find 

oo 

2 f dk, 



x 



J ^k 2 [ln(r( S + A; 2 ))+ 7E ] 



e ik(x 1 -x 2 ) e ik(x!+x 2 ) 
X!>0 x 2 >0 

e -^/s\x-X 2 \ _ e ~Vs(x2+x) 

< Vl 

e —y/s\XQ—Xi\ _ g— ^S(X!+X ) 



2y^ 

Performing the space-integrations, we find 
Z^(x ,x,s) 

CO 

/dk r i r 

— cos(fcxo) - e~^ x ° cos(kx)-e-^ 
2tt L J L 



(48) 



fc 2 [ln(r(5 + fc 2 )) + 7E ] 
X (s + k 2 ) 2 

Note that this is (rescaling k — > ^/sk) 



(49) 



OO 

Z+(x ,x,s) = — / — 
x J 2tt 

— oo 



COS 



(kXy/s) 



-x^s 



cos(kxo^) - e~ x °^ 
k 2 [In (rg(l + fe 2 )) + 7 e 

v^(i + & 2 ) 2 



(50) 



The next step is to invert this Laplace transform which 
is performed in Appendix |B 3| 



variable z = x/y/2t. Note that setting O = 1 — 4e, i. e. 
a\ = —4, the term Z+\x,i) does not depend on x$: 



(z 2 - 2) [\n(2z 2 t) + 7E ] - 2 + I(z) + c(t) 



c(t) = In(t) - 2 7E + 2 



(51) 



where T(z) is defined in Eq. (B53). The first line is ar- 
ranged as to not contribute to the normalization, whereas 
c(t) is independent of z and will not appear in the final 
conditional probability. 7e is Euler's constant. The prob- 
ability distribution, P + to find a non-yet-absorbed 
particle in the interval (x,x + dx) can be computed fol- 
lowing the lines of Eq. ([24]) to order e as 



zf(x,t) + cZ^\x,t) 
Jdx (z^ 0) (x,t) + eZ^ 1) (x,t))) 



oo 

jdxZf(x,t) 





/ 00 r \ \ 



Z ( +\x,t) J dxZ f\ x ,t) j 
/ 



(52) 



Note that the term proportional to c{t) cancels in nor- 
malized objects such as P + (x,t). Therefore, we obtain 

P+(x,t) dx = R ( °\z)dz 



{l + e [(z 2 - 2) ( 7E + ln(2z 2 r)) - 2 • 



(53) 



where R+\z) = zexp(— z 2 j2\ and T(z) is given in 
Eq. ( B53| ). The result in Eq. ( [53] ) still involves both z 
and t. The reason is that for H ^ 1/2 the natural scal- 
ing variable is y = x/(\^2t 1 / 2 + e ) instead of z = x/y/2t, 
as can be seen from Eq. ([2]). To rewrite Eq. (53) in terms 
of y = zt e , we note that 



R i °\z)dz = R { + ) (yt e )t e dy 
= 4%)(l + e 



(o), 



yd y Rf(y) 



i 



hit 



= R^(y){l-e[y 2 -2]\nt} . 
This gives for Eq. ( |53| up to terms of order e 2 
P+(a;,t) dx = Rf(y)dy 



1 



(y 2 -2) ( 7E + ln(2y 2 )) 



(54) 
(55) 



This is the final result announced in equation (|9|, with 
X(y) calculeted in (B53) and below. 



D. Step 4: The Probability P+(x,t) 



Th e final result f or Z ±\x, t) is given in Eqs. (B60) and 
(B65) of Appendix |B 4| expressed in terms of the scaling 



V. COMPARISON TO NUMERICS 

In this Section, we compare our analytical results with 
numerical simulations. More specifically, we consider the 
super- diffusive process with H = |. 
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R + (y) R + (y) 



l r 




FIG. 5: (Color online) R+(y). Analytical result for Brown- 
ian motion R+\y) = ye~ y ^ 2 (solid blue line), and simulation 
data for L = 20000, and H = - (red dots), as well as for 
the fBm with H = | (black diamonds). Histograms are per- 
formed over 4 x 10 5 paths. 



FIG. 6: (Color online) R+(y). Analytical result for Brown- 
ian motion R+\y) = ye~ y ^ (solid blue line), and simulation 
data for L — 20000, and H = | (red dots), as well as for 
the fBm with H = | (black diamonds). Histograms are per- 
formed over 4 x 10 5 paths. 



A. Methodology of simulations 

We aim to sample a fBm processes x(t) at discrete 
times t\ = 1,^2 = 2, . . . , = L. The covariance matrix 
of . . . , Xi, . . . , xl) coincides with the autocorrelation 
function of the original fBm process in Eq. ([6]), setting 
D = l, 



{X{X j 



i 2H +j 2H 



■J 



\2H 



(56) 



The L x L covariance matrix C is symmetric and has 
positive eigenvalues; it is thus possible to find a matrix 
A, positive and symmetric, such that C = A 2 . Matrix A 
is called the square root of C. 

One can simulate paths of a fBm using the stan- 
dard procedure for Gaussian correlated processes: (i) 
Determine A, the square root of C. (ii) Each path 
x = {xi, . . . , Xi, . . . , xl} is given by the matrix multi- 
plication x = Aff. The vector ff = {771, 772, ••• , tjl} is 
a set of L independent Gaussian numbers with unitary 
variance and zero mean. It is easy to check that these 
paths are characterized by the correct covariance matrix 
{56). 



Unfortunately this procedure is time consuming, as for 
step (i) it requires the full diagonalization of C. Better 
results are obtained by making use of the stationarity of 
the increments & = Xi — (we set x\ — £1). Using 
Eq. (56) we can compute (5, the covariance matrix of the 



increments, 



Ci, i+k : = m +k ) = ifc-ir +(k+iy H -2^ , (57) 

where k = 0, . . . , L — i, and Ci+k,i = The matrix 

C is symmetric and positive definite like the matrix C, 
but it also is a Toeplitz matrix. For Toeplitz matrices 
efficient numerical methods allow to avoid the full diag- 
onalization of C. In particular, the Levinson algorithm 



(for a practical implementation of Levinson's algorithm 
see [39] and [40]) is suitable for first passage problems, as 
it recursively generates the increment <^ + i given £1, 
The points of the fBm path are given by xi = Yl]=i £r 
In our simulation we are interested only in positive paths 
(xi > for all i). The Levinson method allows to discard 
negative paths whenever a Xi < is generated, without 
building the full path. 



B. Simulation results 

For each positive path we record the final position xl- 
The histogram of the rescaled variable y := xl/(2L h ) is 
the scaling function R+(y). The results for H = 2/3 and 
the Markovian case H = 1/2 are presented on Figs. [5] and 
[6] For small y the scaling function, R+(y) behaves as a 
power-law, with an exponent (j). For H = 1/2 we expect 
<j) = 1, for H = 2/3 we expect <p — 1/2. Inspired by 
our perturbative calculation we predict that for y — >> 00, 
R+(y) behaves like ~ y 1 e~ v I 2 . In order to facilitate the 
comparison, we define the scaling function 



r+(y) := e* R+(y) . 



(58) 



The numerical data for the scaling funtion r+(y) defined 
in Eq. (58) are shown on Fig. [F] for H = 2/3. They 
clearly show two distinctive power-law behaviors: For 



small y this power law is 



with 



2 ' 



predicted 



by the scaling relation (j) = ^T^- F° r large y a larger 
exponent 7 = 0.7 ± 0.03 is measured. This is consistent 
with the perturbative calculation, which suggests 7 > (j) 
for H > 1/2 and 7 < (j) for H < 1/2. 

A more accurate comparison between the numerical 
data and the perturbation theory is possible. Our per- 
turbative result given in Eq. ([9| is equivalent to r+(y) = 
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FIG. 7: (Color online) Th e numerically determined function 
r+(y), denned in Eq. (58) for the fBm with H = | (black 
diamonds), using L = 20000 and 4 x 10 5 paths. The asymp- 
totic small- y behavior is consistent with <j) = -^p- = \- The 
large-?/ asymptotics (including the amplitude) was taken from 
(59ft, with slope 7 « 1 - 2e = f . 



j/ [1 + eW%) 
we use 



0(e 2 )] . In order to compare to numerics, 



r c + (y) =ye eW ^ +0(e 2 ) . 



(59) 



While the two expressions are equivalent to order e, the 
latter ( 59 ) has the merit to resum the logarithms for small 



and large y into the power-law behavior 



for y — > 
for y — > 00 



where the exponents are the order-e results 
4> e = 1 - 4e , 7e = 1 - 2e . 



(60) 



(61) 



2/3, i.e. e = 1/6, we predict a scaling ~ 



For # 

7 e = |, using (|61|). Note that the curve drawn is exatly 
the asymptotic oehavior of our analytical result (59), us- 
ing ( [To] ), thus also the amplitude and not only the ex- 
ponent are estimated. This can more clearly be seen on 
Fig. |8| where the solid (blue) line represents the theoreti- 
cal order-e prediction, and the dashed line the asymptotic 
behaviors given in Eq. (60). 



Conversely, relation (59) can be used to extract W(y) 
from r_|_(2/), 



FIG. 8: (Color online) Blue solid line: The function r+(y), 



defined in Eq. (59) for H = 



are the predicted asymptotic behaviors, 



The dashed lines 
* y^ (for small y) 
and rsj y 1 (for large y). Superimposed are the simulation data 
shown on Fig. [Tj Note that there is no fitting parameter. The 
deviation at small y is due to the fact that e is rather large, 
so the order-e slope (j) ~ 1 — 4e = | is smaller than the exact 
result 4> = = 1/2+e > wn i cn evaluates to ~. For large y, 
but smaller than 3, the effective cutoff of the simulation, both 
amplitude and slope are correctly predicted. 



W(y) 




logOO 



FIG. 9: (Color online) Blue solid line: The function W(y), 
defined in Eq. |To| . Black diamonds: Estimation of W(y) 
from the numerical data for H = | , using relation ( 62 ) . The 
agreement is quite good for 1 < y < 3. It breaks oown for 
y > 3 due to numerial problems. For y < 1, the deviations 
can be attributed to the large value of e. 



VI. CONCLUSIONS 



^)^m(^) • (62) 

This relation should work the better, the smaller e is. 
Using our numerical results for H = |, we obtain the 
curve presented on Fig. [9j The agreement is quite good 
for 1 < y < 2.5. It breaks down for larger y due to 
numerial problems. For y < 1, the deviations can be 
attributed to the large value of e. 



In this article, we develop a systematic scheme to cal- 
culate the corrections to the universal scaling function 
R+(y) for fractional Brownian motion, in an e = H — \ 
expansion. We compute the full scaling function R+(y) 
to first order in e. In particular we find that R+(y) 
behaves as R+(y) ~ y^ as y — > (near the absorbing 
boundary), while R+(y) ~ y 1 exp(— y 2 /2) as y —> 00 
(far from the boundary), with, at the first order in e, 
(j) = 1 - 4e + 0(e 2 ) and 7 = 1 - 2e + 0(e 2 ). For 
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small e our results confirm the scaling relation found in 
Ref. [29 : R+(y) ~ y* with <j> = 0/H. For fractional 
Brownian motion it is known that = 1 — H, so that 
(j) = (1 — H)/H ~ 1 — 4e+. . .. Far from the boundary, i.e. 
for large y, the leading behavior R+(y) ~ ex.p(—y 2 /2) re- 
covers the Gaussian propagator Q in absence of bound- 
aries; our approach shows that R+{y) has a subleading 
power law prefactor y 1 where 7 is a new (independent) 
exponent. 

Our numerical simulations show that the predictions of 
the asymptotic behavior of R+(y) hold at H = 2/3. In 
particular the two exponents 7 and (j) have been measured 
and shown in Fig. [7] 

Let us stress that few results are known about non- 
Mar kovian processes in presence of boundaries. Pertur- 
bation theory developed in this paper can provide sub- 
stantial new insight here. The method is versatile and 
can in principle be extended to the calculation of other 
quantities such as the propagator for a process confined 
to a finite interval with absorbing boundaries, or alter- 
natively with other, e.g. reflecting boundary conditions. 
Particularly interesting for applications would be the hit- 
ting probability Q(x,L), the probability that a generic 
stochastic process starting at x and evolving in a box 
[0,L] hits the upper boundary at L before hitting the 
lower boundary at [41 j. In the context of polymer 
translocation, the hitting probability is the probability 
that a finite polymer chain will ultimately succeed in 
translocating through a pore. 

In the more general framework of anomalous diffusion, 
the presence of boundaries has been especially studied 
for non- Gaussian processes. For instance, Levy flights 
are Markovian superdiffusive processes whose increments 
obey a Levy stable (symmetric) law of index < \i < 2. 
The Hurst exponent is H = \/\i [42J. By virtue of the 
Sparre Andersen theorem [43 , the persistence exponent 
is — 1/2, independent of \i. The Laplace Transform of 
the scaling function R+ (y) has been computed in [44] for 
a generic value of \i. A scaling analysis of this Laplace 
Transform shows that R+(y) behaves as R+(y) ~ y 1 ^ 2 ^ 
as y — >• (this in in agreement with the scaling relation 
(j) — 0/H), while far from the boundary the Levy-stable 
behavior is recovered. 

An increasing interest is devoted to Gaussian processes 
with self-affine anomalous displacements (x 2 (t)} ~ t 2H 
with < H < 1 [Ml I45U471 [49] . Our current results 
apply only to fractional Brownian motion, i.e. self-affine 
Gaussian processes defined by the autocorrelation func- 
tion (|6|. In particular for fBm it is known that (i) the 
process has stationary increments, (ii) = 1— if, and (hi) 
(j) = 0/H. For all other Gaussian processes with Hurst 
exponent H, (i) the increments are non- stationary, (ii) 
7^ 1 — H and we particularly emphasize that, (iii) no 
scaling relation is known between <\> and (unlike in fBm 
where (j) = 0/H). Among such processes it is possible to 
show that the one, defined by the autocorrelation func- 
tion 



describes the subdiffusive behavoir of a tagged monomer 
in an elastic interface which initially was flat [T6] , For 
this process the persistence exponent is known only to 
first order in e [16], whereas neither the exponents </>, 
nor 7 are known analytically. It would be interesting to 
determine the full scaling function R+ (y) for this process 
within our perturbative framework. 
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Appendix A: The action 

The aim of this Appendix is to determine the ac- 
tion <SW[:r], the first correction to the Brownian action, 
<S(°)[x], in the expansion of S[x] in Eq. ( |29| . As a first 
step we expand the autocorrelation function ([6| around 
H = 1/2, setting D = 1. 

{x(t 1 )x(t 2 )) = G- 1 (t 1 M) (Al) 
= [G^]- 1 (ti,t 2 ) + eK(t 1 ,t 2 ) + 0(e 2 ). 

The first term is the autocorrelation function for H = 

1/2, 

[G^]- 1 (t 1 ,t 2 ) = 2mm(t 1 ,t 2 ), (A2) 

the second term gives the correction at first order in e, 

K(t 1 ,t 2 )=2[t 1 ln(ti)+t 2 ln(t 2 ) - \t x - t 2 \ In \t x - t 2 \] . 

(A3) 



Inverting Eq. (Al) and expanding up to order e one gets 



G = G<°> + eGW + 0(e 2 ) 
G (D = - G (o) KG Q) ; 

where G^\ti,t2) is defined as 



(A4) 
(A5) 



dt , G ( °)(ti,t , )[G (0) ]- 1 (t , ,t 2 ) = 5(h -t 2 ). (A6) 



One can check that the kernel of the Brownian action, 
S^[x], i. e., 



G m (h,t 2 ) = -l5"(t 1 -t 2 ) 



(A7) 



satisfies Eq. (|A6|), namely 



1 



&t' 5"(t 1 -t')[G^]-\t'M) 



(x{t 1 )x{t 2 ))~(ti+t 2 



,2H 



t 2 



\2H 



(63) 



= - / dt'<f"(*i - t')mm{t' ,t 2 ) 
Jo 

= -d\ min(ti,t 2 ) 5(h ~ t 2 ) • 



(A8) 
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It remains to compute the term G^. Integrating by 
parts one has 



G w (t u t 2 ) 

= -\ \ [ At' At" 5' (h - t')5'(t 2 - t")d t ,d t „K(t',t") 
4 Jo Jo 



t r t 



JO 



At' At" S'ih-t'^'fo-t") 

xd t ,d t „(\t' -t"\\n\t' -t"\) , (A9) 



using that the first two terms in (A3) do not contribute 



since they only depend on one of the times. The deriva- 
tive is 



d t ,d t „(\t' -t"\\n\t' -t"\) 
1 



W - t n 



25(t' -t")(l + \n\t' -t"\) . (A10) 



The second term is not well-defined. We decide to in- 
troduce a regularization for coinciding times t = t' — » 
In \t — t'\ = lnr where r > should be thought of as the 
time-discretization of the path-integral. Let us first give 
the final result, before commenting on this approxima- 
tion: 



= -ll!l! dt ' dt '' 6 ' {h - t ' ) w^\ s ' {t2 - t '' ) 

-2(l + lnr)[G (0) ]. (All) 
This yields for the action 



f 1 

dh I dt 2 -x(t 1 )G {1) (t 1 M)x(t 2 ) 



JO 

t r t 



dt 



d tl x(ti)dt 2 x(t 2 ) 



\h-t 2 \ 



-2«S (0) [a;](l + lnT). 



(A12) 



We see that the only possibly ambiguous term, the term 
of order lnr, is proportional to the zeroth-order action 
<S(°)[x], thus equivalent to a change in the diffusion con- 
stant D. Thus its effect is easy to check in the final result, 
when looking at observables in a domain unaffected by 
the boundary. 



Appendix B: Evaluation of Z+\x,t) 
1. Evaluation of Zf (x, t) 



This term is easily evaluated. Indeed, Eq. (36) can be 
recast in the following form 

Z*(x,t) 
= -2(1 + lnr) 



r 1 9 

x hm — — 

xo^o xo da 
-2(1 + lnr) 

x hm — — — 
xo^o xo da 

(1 + lnr)- 



/ 

a=l J x 



x{t)=x 
(0)=x 

a 
Art 
x* 
2t 



V[x]e- aS(0) ^e[x] 



3 . 



- t {x+x y 



*ar - --• <B1) 

In going from the first to the second line we have used 
the expression of the propagator in the Brownian case in 
Eq. (20), introducing the factor of a from the observation 
that the latter appears together with x 2 , and readjusting 
the normalization. 



In terms of the variable z = xj v 2t this gives 

Z$(z,t) = zW(z,t)A(z) (B2) 



where Z^\z,t) = ze~ z ^ 2 / (V2^t) is defined in Q and 
A(z) = (1 + lnr) (z 2 -3) . (B3) 

2. (xo, x, s): The integration over k 



We split Zf(xo,x,s) into two parts 
Z+(x ,x,s) = h(x ,x,s) +I 2 (xq,x,s) 



(B4) 



h{xQ,x,s) 



dk 

xq J 2tt 



s(kx^) 



-x^fs 



I 2 (x ,x,s) = — / 
xo J 



dk 

2^ 



cos(kxoVs) - e~ XoV * 

k 2 ln(l + /c 2 ) 
V^(l + ^ 2 ) 2 

cos(kx^) - e~ x ^~ s 

cos(kx \^) - e~ x °^ 

k 2 [ln(rs) + 7e] 
V^(l + ^ 2 ) 2 



(B5) 



(B6) 



a. ii(xo,x,s) 

The expansion of this term for small xq must be done 
with care; when xq acts as a regulator, one cannot simply 
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expand in it. We claim, and show below that 

h{x ,x,s) = lf(x ,x,s) + /f (x,s) + O(x ) (B7) 
with 

4e~ x ^ f dk 



Xq\^S J 2lT 



[cos(/cxqV^) — l] 



oo 

h B (x,s)=4 J 



dk 
2n 



fc 2 ln(l + fc 2 ) 
X (1 + fc 2 ) 2 

cos(kxy/s) - e~ x ^ 

k 2 ln(l + fc 2 ) 
X (1 + fc 2 ) 2 



(B8) 



(B9) 



In order to prove this, we group the four terms in (B5) 
into two times two terms; the first combination is 



/ — ^ 



x 



2tt 



cos(fcxv^) - e x ^ 



fc 2 ln(l + fc 2 
(1 + fc 2 ) 2 



— +4vs + O(x ) 
x 



CJXJ 

/ 



dfc 

2^ 



cos(fcxv^) — e 



-Xy/S 



k 2 ln(l + fc 2 ) 
(1 + fc 2 ) 2 



= Jf^xo, x, s) + I*(x, s) + O(x ) 
where the divergent contribution is 



(BIO) 



It w (x ,x,s) 



CO 

4 f dk 



. ^sCkxy/s) — 

xq J 2tt L 

fc 2 ln(l + fc 2 ) 
: (1 + fc 2 ) 2 " 



(Bll) 



This expansion in xq is justified since ^— ^ — stands out- 
side the integrand, thus does not act as a regulator. 



The second contribution to (B5) is 

oo 

/ 



_4 f dk 
x J 2tt 



(kXy/s) 



COS 



(kxoVs) 



k 2 ln(l + fc 2 ) 



CO 

2 f dk 



xq J 2tt 

-co 



(1 + fc 2 ) 2 

cos(k(x + xo)V^) + cos(fc(x — xo)v^) 
fc 2 ln(l + fc 2 ) 



2e~ xVs cos(kxoy^) 



(B12) 



(1 + fc 2 ) 2 

Since x ^> xq, we can Taylor-expand cos(fc(x + xo)y/s) 



and cos(fc(x — xo)v^) 7 leading to 



dk 
xo J 2tt 

-co 



cos(kx^) — e Xy ^ cos(kxo\fs) 
k 2 ln(l + fc 2 ) 



(1 + fc 2 ) 2 



O(x ) 



= -If v (x , x, s) + /^(xo, a?, 5) + O(x ) (B13) 
The contributions proportional to lf lv cancel between 



(BIO) and (B13), and we arrive at the decomposition 

(B7T 



We now treat the two contributions to (B7). The first 
contribution I^(xo,x,s) can be evaluated analytically. 
After integration over fc we find a Bessel function, which 
can be expanded in xq as 

Ii(x ,x,s) = 



= -4e 



-y/sX 



ln(xo) + ^ ln ( 5 ) + 7e - 1 



(B14) 
+ O(x ) 



The second contribution i;f (x, s) can be evaluated us- 
ing the relation 

fc 2 



(i + fc 2 ) 2ln(1 + &) = 

We rewrite Ii{x, s) as 
I?(x,s)=4 



" d 




d 




du 


u=l 


du 


u=0- 



" d 


d 




du 


u=l du 


u=0- 



It can be split into two parts, 
Ii(x,s) = fio' 
ha(x,s) 

= -4e~ x ^ 

du 

= - [1 + ln(4)] e~ x ^ 
hb(x,s) 



x,s) 


+ hb(x,s 


) 


' d 


d 




du 


u=l du 


u=0- 



" d 




d 






dii 


u=l 


du 


u=0- 





dk 
2^ 



(l + fc 2 )«+ r 
(B15) 



/^fe e ikxyfs _ e —x^/s 
2^ (l + jfc2)«+i ' 

-co 

(B16) 
(B17) 



CO 

/" dfc 1 

y 2^(i+fc 2 )«+ i 

— CO 

(B18) 

ikxy/s^ _|_ £,2^-(u+l) 



(B19) 

To do the fc-integral in /^(x, 5), it is useful to introduce 
the integral representation 

1 r°° 

(1 + fc 2 )-( n+1 ) = r(i + ^ y dz z u e^ 1+k2)z . (B20) 
This gives 
Jib(x,s) = 4 



" d 


d 






du 


u=l dw 


«=0- 





°° d_fc g ifcx\/i 
-™ 2tt 



dzz"e- (1+fe2 ) z 



, (B21) 
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and performing the Gaussian integral over k yields 
hb(x,s) = 4 ' 



" d 


d 






_d^ 


u=l du 


u=0- 





r(i + u) 



Jo 



/2, 



(B22) 



b. I 2 (x ,x,s) 

l2(xo,x,s) can be calculated using residue calculus. 
We use xo < x to expand the expression, choosing every 
pole in the half-plane in which the corresponding expo- 
nential factor converges. The result is 



7e + ln(rs) 
2^/sxq 



y/s(xQ + x) - 1 
y/s(x - x ) - 1 



e ^/~s(x— x ) 



(B23) 



Expanding for small xq yields 



J 2 (x , x, s) = e~^ sx (2 - y/ix) [ln(rs) + 7e ] + O (x ) 

(B24) 



c. Summary of all terms contributing to Z+(xo,x, s) 
It is useful to re-organize 

Z^(x , X, S) = Ii(x , X, S) + ha(x, s) + Ji 6 (x, s) 

+/ 2 (x,s) + O(x ) (B25) 
as the sum of three contributions: 

(x , x, s) = JoOo, ^> s ) + 5 ) + ^0, + O(x ) . 
The first term depends on xq, 

J (x , 5) = e~ x ^ s [3 - 2 7E + 2 ln(r/2) - 4 ln(x ) 
while the other two terms are 

1 



(B26) 



(B27) 



Ji(x, 5) = 4 



" d 


d 






dii 




u=0- 





r(i + u) 



JO 



dz 



/o 2 V^ 

J 2 (x, s) = -xy/se~ Xy/( * [te + ln(rs 



(B28) 
(B29) 



3. Z+(x,t): The inverse Laplace transform of 

Z+(x,s) 



The inversion of Jq(xq, ^, s) is done by observing that 



Z^(x,s) = lim —Z+(x ,x,s) = e~^ x . (B30) 

^ x ^0 Xq 



1 



This yields 



where Z^\x,t) = 2 ^ 3/2 e "** , and 

£ (>o) = 3 - 2 7E + 2 ln(r/2) -41nx . (B32) 

The inverse Laplace transform of the second term can 
be done directly, 



J (x ,x,t) = zf(x,t)B Q (x Q ) , (B31) 





Ji(x,s) 








= 2 


' d 


d 






1 


_d^ 


u=l d ^ 


u=0- 




.r(i + «) 




x J 

J 


3 v 7r 


1 / 2 e- 





(B33) 



We observe that S (j^ — tj = 5 — z^j z/t and obtain 
2 



Ji(x,t) 

Finally 

Ji(x,t) ■ 



' d 


d 




du 


u=l du 


u=0- 



/^2n«+1/2 

V4t ) 

r(i + «) ■ 

(B34) 





1 


X 








' x 2 




X 


2t 





e 4 * 



V 2t 



-2(7E + ln(- 



■In 2 



Introducing the variable z = x/y2t we have: 

z 



(B35) 



Ji(z,t) = -^e-^^B^z) = Z i f(z,t)B 1 (z) , (B36) 
2irt 



where Z^\z,t) = z exp(-z 2 /2)/ v / 27rt, see Eq. (23), and 
£i0) = (z 2 -2)(te - l + 21nz-ln2) -2 . (B37) 
This completes the Laplace inversion of J\(x,s). 

a. J2(x,s) 

The Laplace inversion of the second term J 2 (x,s) is 
more complicated, and we split it as 

J 2 {x,s) = -[ln(rs) +j E ]xVse- x ^ 

= x(\nr + 7 E )^e"^^ + x^- (ins e~ x ^) 



J2a(x,s) + J 2 b(x,s). 



(B38) 



It is easy to perform the Laplace inversion of the first 
term: 

«x,() = (lnr + 7E )-|=A Ie -5j 

-^■■'^^O-s) ■ (B39) 
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The inverse Laplace transform of the second term can be 
written as 



J 2b (x,t) = x — f(x,t) , 



(B40) 



where 



POO 

/ e- st f(x,t)dt = e-^lns = h(s)g 2 (s) ( B41 ) 
Jo 



~ / \ Ins 



(B42) 
(B43) 



The idea is to inverse-Laplace transform gi(s) and g 2 (s), 
and then to calculate f(x,t) as convolution of and 
#2(£), using (41). These inverses are 



1 fx 2 



2V^ \ 2t 
ln(4t) + 7 e 



01 W = 

The convolution is 

f(x,t)= f gi(t')g 2 (t-t')dt' 
Jo 

rt dt' 



1 e 4t 



(B44) 
(B45) 



(B46) 



r dt' 

Jo 2nt'^ 2 W 



1 e 4t 



4 ln(4[t-t / ])+ 7E 



Using (B40) we have 

.2 



X 2 yt dt / 
, * ) "47T i /o ^ 



2t' 



e 4t- 



^ ln(4[t -£'])+ 7e 



(B47) 

Making a change of variables t' = lit, and using z = 
x/y/2t, this gives 



du 



^z 2 



27TtJ U 5 / 2 Vl^ 

xe"fc [ln(4t) +7e + ln(l - u)] (B48) 
The integral contains two pieces, which we note 
(ln(4t)+ 7E )F 2 (z)+F 3 (z) 



The first piece is 



.2 /»! 



2tt 

The second integral 

^ 2 ^ du 



t 



du 



(B49) 



27T J %fil 2 yj\ - U V ^ 

" (, 2 -i). 



3 e-fe 



(B50) 



(B51) 



is more difficult, but can be performed using Mathemat- 
ica. A convenient substitution a = z 2 (l/u — 1) allows to 
write 



#)=c-t — , (B52) 



where 



X{z) = — L— [ ^L\n(^-)(z 2 +a)(z 2 +a-3)e- 
V2^z 2 Jo V® \z 2 +aJ K 

2 F 2 (l, 1; |, 3; y ) + tt(1 - * 2 )erfi(*/V2) 



6 

-3z 2 + V^e^z + 2 . 
erfi is the imaginary error-function, 



erfi(x) := 



r 



(B53) 



(B54) 



The hypergeometric function 2-F2 (l, 1; §,3; 2 2 /2) can be 
denned by its series expansion 



r 2 00 



n!(2z 2 ) n 
(2n + 4)! ' 



(B55) 



The error-function and the exponential function can be 
combined in another converging series, 



£ 

n=0 



n^2n+l 



(2n- l)(2n + l)n! 



(B56) 



While problems of numerical precision appear for y > 7, 
we can use the asymptotic expansion 

l(z) = l- 7E -ln(2, 2 )+^-^ + A +0( ,-8) (B57) 



At £ = 7, the relative numerical agreement of (B57) and 
dB53b is about 1(T 6 . 



Note that J °° dz ze z / 2 l(z) = 0, thus X(z) does not 
contribute to the normalization. 



b. J 2 (x, t) = J2a(x,t) + J2&0M) 

The sum J 2 (x,t) = J 2a (%,t) + J 2 b(x,t) can be ex- 
pressed using the variable z = x/y/2t as 



J 2 (z,t) = Z^_ 0) (z,t)5 2 (z,t) 



(B58) 



where ^ 0) (z,t) = ze~** / 2 /(y/2Kt) and 



B 2 (*) = (> 2 - 1) ln(4t/r) +I(z). (B59) 
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4. Summary of all terms 



In summary, 

Z$\z,t) = zf(z,t) [a(z) + B (x ) 
—a\ In xq 



B 1 (z)^B 2 (z) 
(B60) 



where Z^\z,t) = ze~ z ^ 2 /(y/2wt) is defined in Eg ([23 



-z 2 /2 



(B61) 
(B62) 
(B63) 
(B64) 



The t erms in q uestion are given in Eqs. (B3), (B32 
( B37| ) and ( |B59| ), and repeated here: 

A(z) = (1 + Inr) (z 2 -3) 
B (x ) = 3 - 2 7E + 2 ln(r/2) - 41nx 

= (z 2 -2)( 7E -l + 21nz-ln2)-2 
B 2 (*) = (^-1) ln(4t/r)+X(z). 
Their sum is 

A(z) + £ + £1(2) + B 2 (z) - ai lnx 
= {(^ 2 - 2) [\n(2z 2 t) + 7E ] - 2} + X(*) 
-(4 + ai)lnx + c(t) 
c(t) = ln(t) + 2 - 2 7E . (B65) 

The result is arranged such that the term in the curly 
brackets, when multiplied by Z+\z,i), integrates to 

zero, as does Z+\z,i)T(z). The propagator Z+\z,i) 
becomes independent of xo if a\ = —4, equivalent to 
= 1 — 4e + 0(e 2 ). As expected, 0o = 0? see Eq. ( [l5| ). 
Since c(t) only contributes to the (time-dependent) 
normalization, it does not enter the scaling function 

R+(y)- 

On the other hand, the only contribution to the nor- 
malization of the propagator Z+(x,t) comes from c(t). 
Since Z+\x,i) integrated over x equals 1, we conclude 
that the survival-probability is 



S(x ,t) = r 
~ t~ 

in agreement with # = 
of our calculations. 



lntj 

e + 0(e 2 ) (B66) 



l + e^2-2 7E 
— i^". This is a non- trivial check 



Appendix C: Scaling arguments 

Consider a process x(t'), starting at x(0) = xo, and 
arriving at x at time t, without having crossed zero, i.e. 
x(t') > for all t' < t. Denote Z+(xo,x,t) its arrival 
probability density at x. Further denote 



S(x Q 



poo 

, t) := / dxZ+(xo,x,t) 
Jo 



(CI) 



the survival probability or the persistence up to time t. 
At late times and fixed xo, for many processes, this sur- 
vival probability decays algebraically 



where is the persistence exponent [3]. Let us now as- 
sume that the process x(t) is self-afflne. This simply 
means that the process is characterized by a single grow- 
ing length scale x ~ t H where H is the Hurst exponent 
of the process. For example, ordinary Brownian motion 
is a self-affne process with H = 1/2. Since the only 
length scale is x ~ t H , the survival probability S(xo,i) 
is a function of only the scaled variable yo = xo/t H , i.e, 
S(x ,t) = G (f^). In order that S(x ,t) ~ t~ e for large 
t and fixed xo, the scaling function G(y), for small 
must behave as 



G(yo) ~ yt , 



where 



H 



(C3) 



We next define p Xo (x,t) as the conditional probability 
density of finding the walker, given that it has not been 
absorbed at any previous time: 

Z + (x ,x,t) _ Z + (x ,x,t) 
J dxZ+(x ,x,t) S{xo,t) 



Note that following Eq. (22), the probability distribution 



of a non-adsorbed particle is for xq — » 

P+OM) =po(a:,t) = lim T ^+^^-. (C5) 

We anticipate the following scaling form for Z+(xo,x,t) 

Z+ (x ,x,t) = ^F^,^). (C6) 

In terms of t he s cale v ariab les y = x/t H and ?/o = %o/t H 
we get from (C4) and (C6) 



F(y,y ) = G(y )p yo (y) 



(C7) 



where p yo (y) is the conditional probability density (C4) 
expressed in terms of the rescaled variables. In the long- 
time limit, 2/0 — ^ and F(?/, ?/o) can be factorized as 



F(y,yo) ~ y e /H Po(y) = yl /H R + (y). 



(C8) 



5(xo, t) ~ t 



(C2) 



Let us now consider the limit ?/ — »• and suppose that 
Po(?/) = R+{y) ~ 2/^- The process is time-reversible in- 
variant, since its increments are stationary, i.e., a path 
from xo to x forward in time plays the same role as a 
path from x to xo backward in time. As a consequence, 
F(y,y ) is a symmetric function, F(y,y ) = F(y ,y). 
Factorization of probabilities for x and xq to zero and 
symmetry thus implies F(y,yo) ~ (yoy)°^ H and it fol- 
lows the proposed scaling relation </> = 0/H. 
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